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SUMARY 


Details  are  given  of  the  expressions  used  to  describe  the 
aerodynamics  and  kinematics  of  the  Sea  King  Mk.50  helicopter  during 
3teady  flight  and  low  rate  manoeuvres  up  to  an  advance  ratio  of 
0.3.  The  aerodynamics/kinematics  formulation  is  a  major  component  of 
the  Sea  King  mathematical  model  developed  by  Aeronautical  Research 
Laboratories  (ARL)  for  flight  simulation  of  this  Anti-Submarine 
Warfare  helicopter. 
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NOTATION 


Moments  of  Inertia  of  helicopter  about  its 
respective  body  axes. 

Effective  cyclic  pitch  control  Inputs  with 
respect  to  shaft  axes  after  allowing  for 
control  phasing  and  blade  lag. 

Effective  cyclic  pitch  control  inputs  with 
respect  to  wind  axes. 

Total  cyclic  pitch  variation  observed  in  wind 
axes  including  pitch  flap  coupling,  i.e. 
orientation  of  ANF  axes  with  respect  to  wind 
axes. 

Blade  tip  loss  factors  for  main  and  tail 
rotors. 

Profile-drag  coefficient  of  the  vertical  tail 
at  an  angle  of  attack  of  90  degrees.  (See 
14,  equation  3). 

Profile-drag  coefficients  of  the  fuselage 
when  moving  in  the  xH,  yH  and  zH  directions, 
respectively.  (See  54,  equations  1,2  and  5 
respectively) . 

Coefficients  of  main  rotor  force  components 
along  the  negative  senses  of  the  0xANfr,  0xH, 
0xs  and  0xw  axes,  respectively. 

(Cfl  -  Main  Rotor  Drag  Force/qs ) . 


NOTATION  (CONT.) 


Coefficients  of  main  rotor  moment  components 
(positive  for  right-handed  sense)  about  the 
0xANF’  0xH  and  0xw  axe3»  respectively. 

( Cj^  ■  L/  QgR ) . 

Coefficients  of  main  rotor  moment  components 
(positive  for  right-handed  sense)  about  the 
0yANF*  0yH  and  0yw  axes  respectively. 

(CM  -  M/qgR) 

normal  force  coefficient  of  horizontal  tail 
plane.  (See  i  4,  equation  6). 

Coefficients  of  main  rotor  moment  components 
(positive  for  right-handed  sense)  about  the 
0zANF*  0zH*  0zS  and  0zw  axes»  respectively. 
(CQ  -  Q/qsB) 

tail  rotor  torque  coefficient. 

(CQt  “  QT/qSTRT> 

Coefficients  of  main  rotor  force  components 
along  the  negative  senses  of  the  0zANF,  0zH, 
0zs  and  OZy  axes,  respectively. 

(CT  •  Main  Rotor  Thrust/qs) 

tail  rotor  thrust  coefficient. 

(Cj^  -  Tail  Rotor  Thrust/qgT) 

Coefficients  of  main  rotor  force  components 
along  the  positive  senses  of  the  0yANF,  0yH, 
Oyg  and  0yw  axes,  respectively. 

(Cy  -  Main  Rotor  Side  Force/ qs) 

Factor  relating  downwash  at  fuselage  to 
induced  flow  through  the  rotor  disk. 
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fgloss 

°RAT 

HROT 

(hrot*o 

*T 

Ktv 

L,  M,  N 

Lpt  ^v*  gth,gt 


Mf 

mth 


Factor  representing  losses  in  transmission 
system  and  accessory  drives. 

Ratio  of  tail  rotor  rotational  speed  to  that 
of  main  rotor. 

Height  of  rotor  above  ground. 

Height  of  rotor  when  aircraft  is  on  ground. 

Total  effective  inertia  of  rotors  and 
transmission  about  main  rotor  shaft. 

Constant  used  in  calculations  of  vertical 
tail  force. 

Total  aerodynamic  moments  about  respective 
0xH*  0yH*  0zH  60(1  y  axes* 

Rolling  moment  contributions  arising  from 
fuselage,  vertical  tail  fin,  horizontal  tail 
and  tail  rotor. 

Pitching  moment  contribution  arising  from 
fuselage. 

Pitching  moment  contribution  arising  from 
horizontal  tall. 


Mjjp  -  Blade  tip  Mach  number. 

Np  -  Yawing  moment  contribution  arising  from 

fuselage  and  tail  plane. 


N 


r 


Yawing  moment  contribution  due  to  yaw  rate 


<3»/3rHEl)) 


NOTATION  (CONT.) 


nt 


Tawing  moment  contribution  arising  from  tail 
rotor . 


OLD  A1S 
OLD  B1S 
OLD  e 

c 

Q 


qCMP 

qeng 

ULOAD 


QMAX 


Q* 


Control  system  cyclic  and  collective 
pitch  outputs. 


Aerodynamic  torque  acting  on  main  rotor  about 

its  07  axis;  positive  retarding. 

^ANF 

Compressibility  torque*  main  rotor. 

Torque  supplied  by  engine. 

Total  load  on  engine  arising  from  main  rotor, 
tail  rotor,  accessories  and  transmission 
losses. 

Engine  torque  at  maximum  power  and  specified 
speed . 

Engine  torque  as  registered  by  aircraft 
instrument. 


QT  -  Aerodynamic  torque  acting  on  tail  rotor 

shaft . 

Compressibility  torque,  tail  rotor. 

Tip  radius  of  main  and  tail  rotor 
respectively. 

(RPM)ggj  -  Specified  rotor  r.p.m.  at  which  maximum 

engine  power  is  delivered. 


^T^CMP 
R  $  Ry 


NOTATION  (CONT.) 


Maximum  engine  power. 

Areas  of  horizontal  and  vertical  tail 
surfaces. 

Coefficient  of  quadratic  terms  in  expressions 
for  fuselage  longitudinal  force.  (See  51*. 
equation  1 ) . 

Equivalent  flat  plate  drag  contribution 
arising  from  undercarriage  and  externally 
carried  weapons  (Torpedoes).  (See  51*, 
equation  1). 


Coefficient 

expression 


of  angle  of  attack  term  in 
for  fuselage  normal  force: 
(See  IH,  equation  5). 


Equivalent  pitching  moment  coefficient 
contributions  arising  Trom  fuselage,  under* 
carriage  and  weapons.  (See  5^,  equation  10). 

Projected  areas  of  fuselage  in  body  axes 
directions. 

Components  of  dunking  sonar  cable  tension  in 
body  axes. 

Helicopter  airspeed. 

non-dimen3ionalized  helicopter  airspeed. 
Resultant  airspeed  at  tail  rotor  hub. 
Non-dimensional ized  airspeed  at  tail  rotor 
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NOTATION  CCONT.) 


Horizontal  wind  speed  relative  to  Earth. 

Resultant  forces  along  body  axes. 

Components  of  aerodynamic  force  on  the 
fuselage  due  to  motion  of  the  helicopter 
along  the  respective  body  axes. 

Effective  side  force  arising  from  tail  rotor 
thrust  after  allowance  for  fin  blockage. 

Side  force  contribution  due  to  vertical  tail 


Normal  force  contribution  of  horizontal  tail 
plane . 

Lift  curve  slopes  of  blade  sections  of  main 
and  tail  rotor  respectively. 


Coning  angles  for  main  and  tail  rotor. 

Inertial  acceleration  components  in 
respective  body  axes. 

Speed  of  sound. 

Lift  curve  slope  of  vertical  tail  fin. 

Longitudinal  and  lateral  cyclic  flapping 
coefficients  of  main  rotor. 

Number  of  blades  in  main  and  tail  rotor. 

Distance  of  horizontal  tail  plane  centre  of 
pressure  from  roll  axis. 
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NOTATION  (CONT.) 


empirically  based  constant  used  in 
expressions  for  fuselage  normal  force. 

Ground  effect  factor  used  in  main  rotor 
aerodynamics. 

Function  of  fuselage  advance  ratio,  used  in 
estimation  of  rotor  downwash  on  horizontal 
tail  plane. 

acceleration  due  to  gravity. 

Respective  heights  of  main  and  tail  rotor 
hubs  above  c.g. 

Height  above  c.g.  of  point  of  application  or 
fuselage  side  force;  used  in  rolling  moment 
estimation. 

Height  of  vertical  tail  centre  of  pressure 
above  roll  axis. 

Incidence  setting  of  horizontal  tail  plane; 
positive  for  leading  edge  up. 

Factor  applied  to  longitudinal  flapping  angle 
(a1 ) . 

Induced  velocity  parameter  associated  with 
main  rotor. 

Induced  velocity  parameter  associated  with 
tail  rotor. 

Coefficients  representing  pitch  change  due  to 
blade  flapping  and  lagging  respectively. 


NOTATION  (CONT.) 


Constant  relating  steady  state  blade  lag  with 
rotor  aerodynamic  torque. 

Distance  of  tail  rotor  from  yawing  axis. 

Distance  of  horizontal  tail  plane  centre  of 
pressure  from  pitch  axis. 

Distance  of  vertical  tail  centre  of  pressure 
from  yaw  axis. 

All-up  mass  of  helicopter. 

Angular  velocity  components. 

Dlmensionalizing  factors  applied  to  force 
coefficients  of  main  and  tail  rotors, 

qg  -  pHR2  (OR)2  and 
qgT  -  pHRT2  (fiTRT)2 

Linear  velocity  components. 

Effective  normal  velocity  at  horizontal  tail 

plane. 

Normal  velocity  component  induced  by  rotor 
downwa3h. 

Position  co-ordinates  of  helicopter  with 
respect  to  eartn  referenced  axis. 

Co-ordinate  of  c.g.  location  forward  of 
fuselage  datum. 


NOTATION  (CONT.) 


Distance  benind  c.g.  of  point  of  application 
of  fuselage  side  force;  used  in  yawing  moment 
estimation. 

Matrices  defined  in  Appendix  2  representing 
axis  transformation  through  angles 
*HE’  6 HE’  *HE  respectively. 

Shaft  angular  velocity  of  main  and  tail 
rotors . 

Fuselage  angle  of  attack. 

Angle  of  attack  between  apparent  wind  vector 
and  no-feathering  planes  of  main  and  tail 
rotors  respectively. 

Airframe  incidence  for  minimum  drag. 

Airframe  incidence  for  zero  lift. 

Incidence  w.r.  to  minimum  drag  incidence, 
0  "  aD 

Effective  fuselage  angle  of  attack. 

Flow  angle  of  attack  for  horizontal  and 
vertical  tail  surfaces. 

Fuselage  angle  of  sideslip,  also,  blade 
flapping  angle. 

Lock  numbers  for  main  and  tail  rotors . 


Mean  profile  drag  coefficient  for  main  and 
tail  rotor  blade  sections  (both  NACA  0012). 


NOTATION  (CONT.) 


63.  a, 


^0*  ^OT 


cr  cit 


Angles  associated  with  respective  flapping 
and  lagging  hinges  which  give  coupled  pitch 
responses  to  blade  flapping  and  lagging 
motions. 

Flap  hinge  radius  as  a  fraction  of  blade  tip 
radius  for  main  and  tail  rotors. 

Blade  root  radius  as  a  fraction  of  blade  tip 
radius  for  main  and  tail  rotors. 

Effective  main  rotor  collective  pitch  at  zero 
hub  radius  (control  system  input  is  OLD  6C). 

Tilt  angle  of  main  rotor  shaft  relative  to 
body  axis,  tilt  normally  negative. 

Tail  rotor  collective  pitch  input  from 
control  system. 

Effective  tail  rotor  collective  pitch. 

Blade  twist  between  hub  and  tip;  negative  for 
usual  case  where  pitch  decreases  towards  tip. 

Control  phasing  angle. 

Non-dimensional  total  flow  averaged  across 
the  no-feathering  planes  of  main  and  tail 
rotor;  positive  in  the  positive  sense  of 
rotor  thrust. 

Coefficient  describing  linear  variation  of 
A  in  longitudinal  direction. 


Advance  ratio  for  main  and  tail  rotors 


'HOV'  vHOVT 


NOTATION  (CONT.) 


Fuselage  advance  ratio,  Uj^/QR 

Mean  induced  velocity  at  N.F.P.  of  main  and 
tail  rotors;  positive  in  sense  opposite  to 
positive  sense  of  rotor  thrust. 

Mean  induced  velocity  of  main  and  tail  rotors 
at  hover,  as  given  by  momentum  theory. 
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Non-dimensional  induced  velocity  at  N.F.P. 
e.g.  v  -  v/vHOy. 


Lagged  value  of  v  for  calculation  of 
transient  inflow  ratio. 


Steady  state  blade  lag  angle. 
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(♦,  e,  f). 


Air  density. 


Solidity  of  main  and  tail  rotors. 


Time  constant  for  calculation  of  v. 


Time  constants  in  engine  model, 


Euler  angles  expressing  aircraft  attitude 
with  respect  to  Earth  axis. 

Blade  azimuth  angle  measured  in  direction  of 
rotation  (zero  at  downwind  position). 


Earth  axis  azimuth  angle  defining  direction 
from  which  horizontal  wind  component  comes. 


Azimuthal  rotation  of  W-axis  from  S-axis. 
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NOTATION  (CONT.) 


tweff 

Effective  azimuthal  angle  between  w-axis  and 

axis  defined  by  cyclic  control  inputs. 

XANF 

- 

Wake  sweep-back  angle  in  ANF  axis  system. 

XFA*  XLAT 

- 

Wake  sweepback  components  in  body  axis  fore- 

aft  and  lateral  directions. 

Abbreviations 

ANF 

- 

axis  of  no  feathering 

NFP 

- 

no  feathering  plane 

TPP 

- 

tip  path  plane 

c.g. 

- 

centre  of  gravity 

r .p.m. 

- 

revolutions  per  minute 

Suffixes 

Vector  components  not  specifically  given  in  NOTATION  are 
identified  by  the  use  of  triple  subscripts. 

(a)  The  first  subscript  indicates  to  which  quantity  the  vector 
component  belongs,  using  the  notation 

H  -  helicopter  centre  of  gravity 
T  -  tail  rotor  hub 
W  -  ambient  wind. 

(b)  Tne  second  indicates  the  datum  to  which  the  vector  is  related, 
where 

E  -  Earth  axis 
W  -  ambient  air. 


1 .  INTRODUCTION 


A  mathematical  model  of  the  Sea  King  helicopter  in  the  ASM  mode 
has  been  developed  by  ARL,  in  response  to  a  Royal  Australian  Navy  (RAN) 
Tas»<  requirement.  General  descriptions  of  this  model  nave  been  reported 
in  Refs.  1  and  2.  A  more  detailed  description  of  one  of  its  major 
components,  namely  tne  aerodynamics/kinematics  model  is  the  subject  of 
this  report.  The  other  major  components  are  the  control  systems  and 
cable/sonar  transducer  models  which  have  been  described  in  References  3-6 
and  7 ,b  respectively.  The  earliest  version  of  the  Sea  King  model  resulted 
essentially  from  the  substitution  of  relevant  Sea  King  data  into  Packer's 
model  for  the  Wessex  helicopter  described  in  Refs.  9— 11-  Because  the 
Wessex  model  was  initially  intended  for  simulation  of  flight  behaviour  in 
near-hover  conditions  and  at  low  advance  ratios,  y  -  0.1  or  so,  certain 
simplifications  had  been  made  in  regard  to  the  induced  flow  and  fuselage 
aerodynamics.  For  these  and  other  reasons  to  be  discussed  later,  the 
initial  Sea  King  model  has  been  modified  by  the  introduction  of  new 
features  described  herein,  which  will  meet  the  requirement  for  adequate 
simulation  up  to  at  least  y  -  0.3. 

Generally,  the  essential  structure  of  the  original  program  of 
Ref.  9  has  been  retained  together  with  its  unique  subscripting  system  for 
identifying  vector  and  attitude  components.  In  addition  to  the 
modifications  previously  mentioned,  several  other  features  have  been  added 
to  give  a  better  representation  of  the  Sea  King,  such  as  main  rotor  pitch 
flap  coupling  and  cyclic  control  mixing.  Tne  incorporation  of  a  simple 
engine  model  was  also  found  to  be  necessary  in  order  to  account  for  large 
excursions  of  rotor  speed  and  torque  observed  in  flight  under  certain 
conditions.  An  effort  ha3  also  been  made  to  model  the  high  transient 
loads  arising  from  'inflow  damping'  when  rapid  collective  pitch  changes 
are  introduced  in  flight. 

This  aerodynamics/kinematics  model  is  a  dynamic  three- 
dimensional  representation  of  body  motion  with  six  degrees  of  freedom. 
Main  and  tail  rotor  flapping  motions  are  calculated  but  the  rotor  is 
considered  qua3i-static,  in  that  the  disk  is  assumed  to  respond 
instantaneously  to  the  motion  of  the  fuselage.  The  inclusion  of  cross¬ 
coupling  terms,  such  as  arise  from  angular  rates,  means  that  the  model  is 
not  limited  to  small  perturbation  studies.  Being  derived  from  the  Wessex 
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model  of  Ref.  9,  the  present  model  is  also  expected  to  be  capable  of 
predicting  low  rate  (<  5  rad/s)  dynamic  response  manoeuvres.  The  validity 
of  the  model  is  presently  being  assessed  by  comparing  model  predictions 
with  results  from  A.R.L.  flight  trials  of  the  Sea  King  MkbO  as  reported  in 
Reference  [12]. 

Section  2  of  this  report  gives  a  general  description  of  the  Sea 
King  Math  Model.  The  model  structure  allows  convenient  subdivision  into 
self  contained  topics  which  are  commented  on  in  tne  following  sections. 
Section  3  describes  the  main  and  tail  rotor  aerodynamics,  while  the 
fuselage  aerodynamics  are  given  in  Section  4.  The  equations  of  motion  of 
a  rigid  body  witn  rotating  suo  systems  are  detailed  in  Section  5  and  the 
kinematic  equations  are  stated  in  Section  6.  A  description  of  the  simple 
engine  model  used  is  given  in  Section  7.  followed  by  a  discussion  of  the 
control  inputs  in  Section  8.  Throughout  thi3  report,  discussion  is 
limited  to  the  assumptions,  approximations  and  physical  aspects  related  to 
the  formulation  of  expressions  used  in  the  coding.  Detailed  derivations 
of  expressions  wnich  may  be  found  in  standard  tests  and  references  are  not 
given. 
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2.  GENERAL  DESCRIPTION  OF  MODEL 

The  aerodynamics/ kinematics  model  constitutes  one  of  3  major 
elements  of  the  complete  Sea  King  mathematical  model,  the  others  being  the 
Systems  and  Cable  models  respectively.  The  computer  program  for  the 
complete  model  i3  written  in  CSMP-10  (ARL)  [13  &  14].  This  is  a  block 
oriented  simulation  language  in  which  the  model  is  expressed  in  terras  of 
linked  modules  or  blocks  which  represent  a  particular  function  or 
operation.  The  language  also  incorporates  "used-def ined"  blocks  written 
as  Fortran  subroutines  which  allow  complex  algebraic  expressions  to  be 
handled  conveniently.  The  Aerodynamics,  Systems  and  Cable  models  are  such 
"uaer-def ined"  subroutines. 

A  schematic  diagram  of  tne  Aerodynamics/Kinematics  model  is 
shown  in  Fig,  1.  Motion  of  the  aircraft  is  obtained  by  applying  the 
equations  of  motion  after  all  force  and  moments  contributions  at  the 
aircraft  c.g.  have  been  determined.  These  components  arise  from  the  main 
rotor,  tail  rotor  and  fuseiage/empennage  aerodynamics  which  may  include 
rotor  downwash  interference  effects.  Axis  transformations  are  required  30 


(3) 


that  all  contributions  are  converted  to  a  common  axis  system  before 
summation.  For  instance,  main  and  tail  rotor  quantities  are  initially 
derived  in  the  usual  axes  of  no-feathering  (ANF)  centred  at  the  respective 
rotor  hubs.  Details  of  the  various  axis  systems  used  in  this  program  may 
be  found  in  Appendix  1. 

Aircraft  accelerations  obtained  from  the  equations  of  motion  are 
* 

integrated  to  give  both  linear  and  angular  velocity  components  in  the 
body  axes  system.  The  orientation  of  this  system  with  respect  to  Earth 
based  axes  is  expressed  in  terms  of  Euler  angles.  These  angles  are 
obtained  from  integration  of  the  Euler  rates  given  by  the  solution  of  the 
kinematic  equations.  Transformations  employing  Euler  angles  are  used  to:- 

(a)  derive  aircraft  velocity  components  in  Earth  axe3  whence  by 
integration  the  position  co-ordinates  are  obtained; 

(b)  express  Earth-referenced  ambient  wind  conditions  in  body- 
referenced  wind- velocity  components.  Knowing  the  aircraft 
velocity,  relative  airspeed  components  may  be  deduced. 

Calculation  of  the  required  aerodynamic  quantities  associated 
with  the  main  and  tail  rotors  requires  a  knowledge  of  the  inflow  and 
advance  ratios  for  each  rotor.  These  are  derived  from  the  local  airspeed 
vector  and  its  angle  of  attack  with  respect  to  the  ANF  for  each  rotor. 
Thus,  the  previously  determined  body  airspeed  components  must  be 
transformed  into  components  appropriate  to  each  ANF  axis  system  via  the 
wind  shaft  axes.  The  angular  orientation  of  the  ANF  to  these  axes  is 
determined  by  cyclic  pitch  control  inputs  together  with  pitch 
contributions  arising  from  pitch-flap  coupling  as  explained  in 
Section  8.  In  the  case  of  the  tail  rotor  no  cyclic  control  inputs  can  be 
applied  so  that  all  cyclic  pitch  changes  stem  from  flapping.  Cyclic 
control  inputs  to  the  main  rotor  must  also  be  referred  from  body  to  wind 
shaft  axis  after  due  allowance  for  control  phasing.  (See  Section  8). 


Integrator  blocks  in  CSMP-10  (ARL.)  normally  use  a  second-order  Runge- 
Kutta  method.  An  Adams  predictor-corrector  integration  method  is  also 
available  [1A], 


(4) 


This  latter  feature  is  introduced  in  the  swash  plate/rotor  head  geometry 
in  order  to  give  a  more  nearly  longitudinal  tip  path  plane  tilt  in 
response  to  longitudinal  stick  displacements. 

The  mean  induced  flow  component  of  the  total  inflow  through  the 
rotor  is  based  on  momentum  theory  together  with  the  Glauert  hypothesis  for 
forward  flight  (Ref.  9,  page  21).  It  is  also  assumed  that  there  is  a 
linear  variation  of  induced  flow  in  the  streamwise  direction,  the 
magnitude  of  which  depends  on  the  wake  sweepback  angle. 

Equations  for  the  forces  and  moments  developed  by  flow  over  the 
fuselage  and  tail  surfaces  are  generally  based  on  those  used  in  the  Sea 
King  Flight  Simulator  [15]  but  rewritten  in  non-dimensional  form.  A 
downwash-modifylng  function  dependent  on  wake  angle,  similar  to  that 
employed  in  Ref.  15,  i3  used  to  adjust  pitching  moments  as  a  function  of 
speed,  so  that  agreement  with  flight  values  of  pitch  attitude  is  improved. 

By  its  nature,  the  simplified  engine  model  does  not  require 
incorporation  in  the  aerodynamics/kinematics  "user  defined  subroutine"  anc1 
is  written  directly  in  the  simulation  language  CSMP-IO(ARL). 

3.  ROTOR  AERODYNAMICS 

3.1  Main  Rotor  Aerodynamics 

Using  blade  element  theory,  expressions  for  hub  forces  and 
moments  in  ANF  axes  may  be  obtained  by  integration  of  the  blade  element 
forces  with  respect  to  radius  and  azimuth,  as  given  in  Ref.  9  and  standarc 
texts.  The  assumptions  of  Ref.  9  have  been  retained  and  are  - 

(a)  Blades  are  infinitely  stiff  in  torsion  and  bending  and  hence 
second  order  flapping  can  be  justifiably  ignored. 

(b)  The  effects  of  reverse  flow,  blade  stall  and  compressibility 
effects  are  not  accounted  for  in  the  blade  element  analysis.  At 
the  highest  aircraft  speed  of  120  kn  (y  ■*  0.3)  these  effects  are 
likely  to  become  evident.  (An  empirically  based  compressibility 
correction  has  been  included  in  the  expressions  for  rotor  torque 
given  in  §7). 
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(c)  Small  angle  approximations  are  used  in  respect  of  inflow,  blade 
incidence  and  flapping  motion  angles. 

Final  expressions  for  the  rotor  forces  and  moments,  which  are 
given  later  in  this  section,  employ  non-dimensional  flow  velocities, 
X  and  p,  normal  and  tangential  to  the  no-feathering  plane  (NFP) 
respectively.  Relations  used  to  derive  these  quantities  in  the  model 
program  are  now  given  - 

1 .  Total  airspeed  vector , 


/  (u2  ♦  w2  ) 

y  ^  HWW  HrfW' 


2.  Angle  of  incidence  of  flow  to  NFP, 

w. 


ay  *  artan 


3.  Advance  ratio. 
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p  -  V  cosaw/flR 


4.  Inflow  ratio  (steady  state), 
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t>.  Transient  inflow  ratio  (for  improved  modelling  of  rapid  collective 
pitch  changes), 


where 


X  »  p  tanaw  -  vlag/RR 
;lag  1 
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in  Laplace  form. 


b.  Mean  induced  velocity  v  is  now  calculated;  the  first  step  is  to 
calculate  the  momentum  theory  value  at  hover,  v„n„  , 
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calculate  the  momentum  theory  value  at  hover,  v. 


vHOV  ^  (C^/2) 


Next,  the  non-dimensional  airspeed  vector, 


V  -  V/v. 
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Now,  the  implicit  relation  for  non-dimensional  induced  velocity  is 
solved. 
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v  ■  (v  +  V  -  2Vv  sina  )  2 
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from  whence, 


V  “  ki  V  vHOV/f(z) 


where  kj  is  an  empirically  based  induced  velocity  parameter. 

The  ground  effect  factor  f(z)  i3  calculated  from  the  expression  of 
Reference  16,  namely 
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where  HROT,  the  rotor  height  above  ground  level,  is  given  by. 
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where  -z^EE  represents  the  height  of  the  aircraft  above  the  ground  and 
(HroT)0  is  the  height  of  the  rotor  when  the  aircraft  is  on  the  ground. 

7.  A  longitudinal  variation  of  induced  flow  (Glauert-type)  is  assumed 
using  the  parameter  A1  as  for  Ref.  9  where, 

A1  “  XANF  v/RR 
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Blade  tip  loss  factor* 
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with  k_  being  an  empirically  derived  factor  allowing  for  the 
a1 

discrepancy  between  flight  data  and  theory. 
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10.  Thrust  coefficient. 
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where  B,  approximates  the  loss  in  rotor  thrust  arising  from  reduced 
lift  in  the  blade  tip  region.  Similar  allowances  are  made  to  the 
induced  drag  components  of  the  rotor  in-plane  forces  given  in  the  next 
section.  Expressions  given  later  for  rotor  hub  moments  are  modified 
in  a  similar  manner  using  the  factor  . 


1.  Kotor  Force  coefficients, 
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12.  Resolution  of  hub  force  coefficients  (ANF  axes)  back  to  body  axe3  for 
substitution  in  the  equations  of  motion.  Small  angle  approximations 
are  used  where  appropriate. 

(a)  From  ANF  to  wind  axes  via  the  cyclic  blade  angles,  A1W>  B1W 
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(b)  From  wind 

to 
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azimuthal  angle. 
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c)  From  shaft  to  helicopter  body  axes  via  the  shaft  angle,  0 


13.  Main  Rotor  torque  coefficient  in  ANF  axes  - 
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14.  Main  Rotor  hub  moment  coefficients  in  wind  axes  - 
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1b.  Hub  moment  and  torque  coefficients  must  be  resolved  Dack  to  helicopter 

body  axis  in  a  similar  manner  to  the  force  coefficients.  Because 

torque  is  the  dominant  moment,  and  since  A,  and  B  are  small,  the 

I  W  1W 

torque  coefficient  in  wind  axis  is  assumed  to  be  equal  to  that  in  ANF 
axis  so  that 


With  nub  moments  C.  and  CM  already  in  wind  axes  we  now  resolve 

LW  MW 

successively  from  wind  to  body  axes  through  the  azimutnal  angle, 
and  shaft  angle,  0SH,  to  give 
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3.2  Tail  Rotor 

The  main  outputs  required  from  the  tail  rotor  model  are  thrust 
and  torque.  The  former  provides  a  fuselage  yawing  moment  to  oppose  the 
main  rotor  torque  reaction  and  also  a  side  force  which  must  be  countered 
by  a  lateral  tilt  of  the  main  rotor.  The  tail  rotor  torque  contributes  to 
the  overall  power  demand  on  the  propulsion  system  and  its  reaction  also 
gives  rise  to  a  small  fuselage  pitching  moment.  Other  tail  rotor  forces 
and  hub  moments  are  not  considered  significant. 

The  tail  rotor  differs  from  the  main  rotor  in  the  following 

respects, 

(a)  Blades  are  untwisted. 

(b)  There  i3  no  swash  plate  and  hence  the  control  axis  coincides  with 
the  3haft  axis. 

(c)  The  flap  hinge  angle,  5^  3  4 b  deg,  results  in  a  blade  pitch 
decrease  of  1  deg  for  each  degree  of  positive  flapping. 

This  last  feature  means  that  cyclic  flapping  arising  from 
aircraft  motion  will  produce  corresponding  pitch  variation  with  respect  to 
the  control  axis  which  is  therefore,  by  definition,  not  the  ANF.  Thus, 
the  orientation  of  the  ANF  with  respect  to  the  control  or  shaft  axis  is 
taKen  into  account  when  calculating  the  airflow  angle  of  attack  with 


respect  to  the  ANF  system  (equation  6  below,  also  see  Ref.  9).  Another 
effect  caused  by  a  large  6^  is  tne  change  in  effective  collective  pitch 
with  coning,  (equation  7),  which  i3  much  greater  than  that  for  the  main 
rotor . 

While  thrust  and  torque  expressions  relate  to  the  ANF  system,  it 
13  assumed  that  because  of  the  small  angles  involvea,  the  values  of  thrust 
and  torque  pertain  to  the  hub/shaft  axis  system. 

Calculation  is  essentially  the  same  as  for  the  main  rotor  but 
with  new  axes  systems  centred  at  the  hub  being  defined.  Notation  changes 
involve  the  additional  suffix  "T"  where  required.  Steps  involved  in  the 
calculation  are  summarised  below: 

1 .  Airspeed  components  at  the  hub  arising  from  aircraft  translational  and 
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motion  in 
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axes.  (Downwash  effects  from  the  main 
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2.  Component  of  incident  flow  in  (XOY)T  plane, 
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0.  Airspeed  magnitude, 
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4 .  Advance  ratio, 


VT  cosaT 


5.  Inflow  ratio 


8.  Coning  angle, 

a0T  0T  0TE 

9.  Induced  flow  calculation  (as  for  main  rotor,  ground  effect  ommitted), 
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where  (k^j  is  the  induced  velocity  parameter  associated  with  the  tail 
rotor . 

10.  Blade  tip  loss  factor, 
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Allowance  for  tip  losses  is  included  in  the  following  expressions  for 
tnrust  and  torque. 

11.  Thrust  coefficient. 
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12.  Torque  coefficient. 
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4.  FUSELAGE  AERODYNAMICS 

Aerodynamic  forces  and  moments  on  the  fuselage  and  tail  surfaces 
only  become  apparent  at  high  speeds.  As  the  tail  surfaces  are  positioned 
so  as  to  exert  controlling  moments,  it  is  important  to  model  these  in 
order  to  predict  trimmed  attitudes  which  agree  with  flight  data.  This 
task  is  rendered  difficult  oy  the  fact  that  the  empennage  lies  in  a  region 
of  complex  flow  interaction.  The  main  rotor  downwash  is  non-uniform  over 
the  fuselage,  and  in  translational  flight  the  wake  sweepback  angle 
varies.  In  addition,  the  tall  surfaces  may  be  immersed  in  the  wake  from 
the  forward  fuselage,  including  the  main  rotor  ma3t  and  head. 


(15) 
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Expressions  used  in  the  present  model  for  fuselage  forces  and 
moments  are  mainly  based  on  Westland's  Sea  King  Simulator  [15J,  while  some 
equations  are  retained  from  the  Wessex  program  [yj.  The  former 
expressions  have  been  modified  so  as  to  be  in  dimensionless  coefficient 
form.  These  relationships  are  based  on  limited  wind-tunnel  data  for  a 
fuselage  with  no  rotor.  To  account  for  varying  downwash  influence  in  the 
tail  region,  empirically  based  functions  of  the  wake  sweepback  angle  are 
incorporated  to  give  better  agreement  with  flight  data. 

The  equations  for  aerodynamic  forces  and  moments  acting  on  the 
fuselage  and  empennage,  in  body  axes,  are  stated  below. 

1.  Fuselage  longitudinal  force: 


V  “  2  P  UHWHiUHWH(  CDXF  SX  *  (SCD)oMaM  +  (SCD)u/c+  (S(Vwep 


where  au  -  - -  -  an,  represents  the  incidence  deviation  of  the 

M  UHWH  ° 

fuselage  from  the  minimum  drag  value,  aQ  and  is  limited  within  the 


range  ±0.4  radian. 


iiA  ■ 

!>:>/■ 


«v 

V  v  ■ 

\4  . 

-*  V.  ; 

v.v.v'. 

v.v.vCv 

at” 

.v.y.v.; 

/  v.v.v; 
v.‘< 

•  :-v-iXv 
-  ■  --V 


2.  Fuselage  lateral  force: 
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with  X.  at  beinS  the  lateral  component  of  the  wake  sweepback  angle. 

The  tail  fin  lateral  force  expression  is  made  up  of  two  components; 
tne  first  component  accounts  for  the  unstalled  case  whilst  the  second 
component  applies  to  a  stalled  fin. 
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Lateral  force  due  to  tail-rotor  thrust  is  given  by: 
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and  a  thrust  reduction  factor  of  0.y2b  has  been  applied  to  account  for 
fin  interference  to  the  tail  rotor  (Ref.  1b). 


The  fuselage  normal  force  is  given  by. 
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where  xPa 


,the  fore-aft  wake  sweepoack  angle  is  given  by 


XFA  “  XANF  +  B1VT  C0S^WH  "  0SH 


and  a7  ,  the  effective  fuselage  angle  of  attack  is  given  by 
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Z  UHWH  L  **  UHWH 


and  fjj  is  an  empirical  constant, 


Normal  force  on  the  horizontal  tailplane  is  given  by 
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with  a  ,  the  flow  incidence  at  the  horizontal  tail-plane,  given  by 
TH 


aTH  “  1  +  artan 


and  i  is  the  horizontal  tail-plane  incidence  setting. 


The  effective  tail-plane  life  curve  slope, 


is  given  by  Fig.  2 


while  the  effective  tail-plane  normal  velocity,  w eff,  is  obtained 


Weff  =  WHWH  +  qHEH  ^TH  *  Wi  (  C03°LDB1S  "  °SH)f(px) 

wnere  OLDB1s  is  the  cyclic  pitch  control  input  (see  8)  and  the  normal 
component  of  the  main  rotor  downwash,  w^  is  given  oy 
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EKj^  is  an  empirical  factor  representing  the  magnitude  of  the  downwash 
in  the  region  of  tne  tail  plane 

and  f(ux)  is  given  by  Fig.  3.  where 
"  UHWH/QR 

The  function,  f(px)  »  is  empirically  derived  and  represents  the 
varying  downwash  influence  on  the  tail-plane  as  a  function  of  forward 
speed. 

With  the  exception  of  the  fuselage  pitching  moment,  the  fuselage  and 
empennage  moments  are  expressed  as  products  of  a  force  and  it3 
estimated  moment  arm  from  the  centre  of  gravity.  The  moments  are 
stated  below. 

7.  Tail  rotor  roiling  moment: 

i.-j*  *  Yfh-j. 

8.  Fuselage  rolling  moment: 

LF  '  *F  hF 

y.  Tail  fin  rolling  moment: 

LTV  “  YTV  hTV 

10.  Horizontal  tail-plane  rolling  moment: 

LTH  “  ZTH  dTH 
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11.  Fuselage  pitching  moment  is  given  by. 


Mp  -  -  p  "%WH"vMwHWH~vl  (SCVbOD+  UHWh!UHWhI  { (SCM}u/c  +  (SCM)wep} 


12.  Pitching  moment  due  to  norizontal  tail-plane  is  given  by 


MTH  “  ZTH  ^TH 


13.  Tail  Rotor  yawing  moments 


nt  -  -  yt*t 


111.  The  combined  fuselage  and  tail  fin  yawing  moment  is  given  by: 


NF  ”  -YF  Xyfm  +  Nr  PHEH 


where  the  moment  arm  Xyf-m  varies  with  sideslip  angle  B  ,  as  3hown  in 
Fig.  H,  and  Np,  the  stability  derivative,  (  -  <*N/drHEH  )  as  given 
in  Fig.  5,  is  a  function  of  both  forward  speed  and  the  yawing  sense. 


Pitch  and  yaw  damping  of  the  fuselage  and  empennage  are  included  in 
expressions  12  and  1H  above,  while  other  terms,  such  as  roll  damping, 
are  assumed  to  be  negligible. 


5.  EQUATIONS  OF  MOTION 

In  order  to  derive  the  helicopter  motion  by  use  of  the  standard 
equations  of  motion  for  a  rigid  body  in  helicopter  body  axes,  all  the 
forces  and  moments  derived  in  the  previous  sections  must  be  combined  to 
give  the  total  aerodynamic  forces  (X,Y,Z)  and  moments  (L,M,N). 
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Thus  for  the  forces  we  have  the  following  equations  - 

X  "  XF  "  CHh  %  +  TXH 
Y  -  YF  *  XT  +  XTV  CXH  qS  +  TYH 
Z  “  ZF  +  ZTH  CTh  QS  *  TZH 

where  Tvu,  Tvu,  T,u  are  respective  components  of  the  sonar  cable  tension 
XH  YH  ZH  2  2 

(derived  in  cable  package  [8]),  and  qs  -  p  IIR  (OR)  is  a 

dimensionalizing  factor. 

Similarly  for  the  moments  we  have, 


L  “  qS  (CYH  hr  *  CLH  *  LT  +  LF  +  LTH  LTV 

M  “  qS  (CHh  hr  “  CTh  XCH  +  °Mh  R)  '  CQt  qST  RT  +  MF  +  MTH 

N  -  -CY  qg  xCH  +  Nt  +  Np  ♦  Torque 
H 

where  hp  and  Xq^  are  appropriate  aircraft  dimensions  related  to  respective 
moment  arms  and 

qST  -  pJiRT2  (QtRt)2  for  the  tail  rotor 


Torque,  Q  »  total  aerodynamic  torque  of  the  main  rotor 


qsR 


CMP 


where  QrMP 


is  the  compressibility  torque  of  the  main  rotor  (see 


The  standard  equations  of  motion  are 
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,  N  ,  B^A  .  ♦  21 

HEH  =  C  "  1  C  ‘  PHEH  QHEH  +  B  C 

where  IT  -  the  total  main  rotor  moment  of  inertia  about  the  shaft, 

and  IHUB  -  tne  moment  of  inertia  of  the  main  rotor  about  the 

shaft,  inboard  of  the  flapping  hinges 


The  last  term  in  the  p  and  q  equations  above  accounts  for 
hub  gyroscopic  effects,  while  the  last  term  in  the  r  equation  accounts 
for  rotational  acceleration  of  the  main  rotor  system  (assuming  a  rigid  lag 
hinge;.  Engine  gyroscopic  terms  can  also  be  included  in  the  above 
equations  but  are  generally  negligible. 

The  linear  velocities  (u,v,w)HEH  and  angular  rates  (p,q,r)HEH 
are  found  by  integrating  tne  above  equations. 

inertial  accelerations  are  given  by  the  standard  equations: 
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where 


‘wee 

"VWE 

0o3,wee 

WEE 

"VWE 

Sln’«EE 

and  wWEE  is  given  along  with  VWE  and  H'WEE  as  input  data.  Thus,  the 
required  airspeed  components,  in  body  axes,  are 


u  * 

HWH 

UHEH 

UWEH 

v  * 

V 

v 

HWH 

HEH 

WEH 

w.  .  * 

w  - 

w 

HWH 

HEH 

wEH 

The  inverse  matrix  conversion  is  used  to  find  the  aircraft 
trajectory  by  resoiving  the  inertial  velocities,  obtained  from  the 
equation  of  motion,  from  body  axes  back  to  earth  axes,  thus 
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-  J 

-1 
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V 
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HEE 
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<L 

HE 
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w 

where  the  inverse  matrices  are  given  in  Appendix  2. 

The  position  co-ordinates  (x,y,z)^EE  are  then  obtained  by 
integration. 

The  main  rotor  aerodynamics  (see  $  3.1)  requires  airspeed 
components  and  pitch  and  roll  rates  to  be  specified  in  shaft/wind  axis. 
A3  a  first  step  the  body-referenced  airspeeds  are  resolved  into  shaft  axes 
as  given  below, 

UHWS 
VHWS  ‘ 

WHWS  " 


UHWH  "  WHWH  0SH 
VHWH 

WHWH  +  UHWH  6SH 


The  required  angular  rates  in  wind  axes  are  similarly  derived  to 
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where  small  angle  approximations  are  used  for  the  shaft  tilt  angle, 

Wind  axes  are  now  defined  by  an  azimuthal  rotation  4^  of  the 
shaft  axes  about  the  QZ  axis,  in  this  case  the  shaft  axis,  so  that  the 
incident  wind  vector  is  contained  in  the  XOZ  plane  of  the  wind  axes.  We 
have 


*^~^UHWS  +  V  HWS^ 


(by  definition) 


and  the  azimuthal  angle  is  given  by 


* WH  *  artan  lVHHS/uHWS> 
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Tne  final  quantities  provided  by  the  kinematics  section  are  the 
linear  velocity  components  of  the  airflow  relative  to  the  tail  rotor  hub 
snaft  axes  system  (see  Appendix  1).  These  airspeeds  are  needed  for 
calculation  of  tne  tail  rotor  aerodynamics.  We  have, 
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7.  ENGINE  MODEL 

The  main  requirement  of  the  engine  model  is  to  simulate  the 
variations  in  torque  and  rotor  speed  (r.r.p.ro)  observed  in  flight 
testing.  It  is  of  particular  importance  in  dynamic  response  test3 
involving  collective  and  pedal  inputs  and  of  lesser  importance  for  cyclic 
inputs.  The  collective  inputs  produce  immediate  torque  response  and 
subsequent  fuel  flow  management  to  maintain  the  set  rotor  revolutions  per 
minute  (r.r.p.m).  The  pedal  inputs  produce  yawing  rates  which  are  sensed 
a3  r.r.p.m.  changes  by  the  aircraft-based  tachometer.  Control  inputs  to 
the  cyclic  channels  produce  pitch  and  roll  rates  which  modify  rotor 
loadings.  However,  the  ensuing  rotor  speed  variations  are  much  less  than 
for  the  collective  and  pedal  cases. 

Detailed  modelling  of  fuel  management  and  speed  control  systems 
may  be  quite  complicated  [17]  when  the  stiffness,  inertia  and  damping  of 
all  the  transmission  components  is  allowed  for.  The  present  simplified 
model,  shown  in  block  diagram  Fig.  7,  contains  some  features  of  Padfield’s 
simplified  model  [18]  notably  the  lead-lag  representation  of  the  engine 
response.  However,  some  difficulties  have  been  experienced  in  choosing 
appropriate  values  for  the  time  constants  which  have  been  estimated  by 
trial  and  error. 

The  diagram  is  largely  self-explanatory,  the  following  equations 
defining  quantities  either  input  to  the  model  or  certain  outputs  restated. 

Total  torque  load  on  the  engine,  i.e.  main  rotor,  tail  rotor, 
accessories,  etc. 


Q  a  (u  +  G  y  j  FQ 
WLOAD  RAT  LOSS 


where  Q  =  C  R  Qs  -  QCMp 

O 


and  Q«C  Rq  ♦  (0  ) 

WT  UT  T  qST  vvt'cmp 


with  G„..p  »  FG.  =  factor  representing  losses  to  transmission  and 

nAT  1  LUoo 

accessories . 
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The  torque  contribution  arising  from  compressibility  effects  is 
given  by  the  following  empirically  based  expression  - 

Hover  case  (p  <  0.035) 

Qcmp  -  0.25  AM2  qs  R  o  (1  -  p/0.035) 

where 

AM  -  (M  p  -  0.78)  ♦  1 .25  CT/0 

and 

MTrD  =  Mach  number  at  blade  tip  »  OR/a 
Hr  3 

Forward  flight  case  (p  2  0.035) 

QCMP  =  1/8  (0-2  AM*  +  0,0085  AM)  Qs  R  o  O  ♦  4.7  p2) 

where 

AM  -  (MTIp  -  0.82)  ♦  2.H  CT/o 

Here  MTIp  refers  to  the  tip  Mach  number  of  the  advancing  blade, 
-  (V  ♦  flR)/as 

Similar  expressions  are  used  in  the  tail  rotor  torque  calculation. 

Engine  output  torque  is  redefined  for  comparison  with  flight  data 
(i.e.  aircraft  instrument)  by  the  following 

QJ  ”  (QENG/UMAX)  X  111,0 

where  the  torque  at  maximum  engine  power  and  specified  speed  is  given  by, 


•&V: 

v  “.  V  • 


(SHPmax  X  550)  /  (RPHset  X  g  ) 


follows. 


The  engine  response  time  constants  x^  and  are  specified  as 


x 2  -  0.008  -  0.00045  Q% 


t3  -  0.06  -  0.005  Qj 


Rotor  speed  (r.r.p.m.)  as  sensed  by  the  aircraft  instrument,  is 
tne  sum  of  the  rotor  angular  velocity  and  fuselage  yawing  rate. 

The  reference  r.p.m.  is  dependent  on  the  collective  stick 
setting,  thus  providing  an  approximate  droop  cancelling  effect.  No 
derivative  term  is  presently  used  to  model  the  anticipator  action. 

A  first  order  lag  (time  constant  r1  )  represents  tne  fuel  flow 
response  to  r.p.m.  errors,  whilst,  as  stated  previously,  the  engine 
response  is  represented  by  a  lead-lag  term. 

The  total  fuel  flow  output  is  an  approximation  to  values  given 
in  the  Sea  King  Operating  Data  Manual  (0DM)  [!<#]. 

8.  CONTROL  INPUTS 

As  with  mo3t  helicopters,  the  Sea  King  incorporates  a  measure  of 

control  phasing,  mainly  to  reduce  the  degree  of  coupling  between  the 

aircraft  pitch  and  roll  response  to  pure  longitudinal  or  lateral  cyclic 

control  inputs.  This  phasing  is  achieved  by  the  geometric  arrangement  of 

the  pitch  link  attachments  to  tne  swash  plate  and  blade  cuffs.  For  the 

Sea  King,  a  value  of  phase  angle,  <  -  17.5  deg  is  used,  wnich  means  tnat 

for  a  fore-aft  jack  movement  the  blade  pitch  changes  are  greatest  when  the 

blade  is  at  an  azimutnal  location  of  »  V0  +  17.5  »  107.5  deg. 

In  the  present  model ,  an  allowance  has  been  made  for  the  steady  state 

2 

blade  lag  angle,  £,  (  »  k^Q/p  )  which  reduces  the  effective  phase  angle. 
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The  systems  package  of  the  Sea  King  model  [3]  generates  cyclic 
blade  angles  A1S  and  B1S  which  are  input  to  the  Aerodynamics  package 
presently  described.  This  notation  has  been  retained  from  the  Wessex 
model  where  no  allowance  for  control  phasing  was  made  and  A1S,  81S 

therefore  truly  represented  the  cyclic  pitch  variation  with  respect  to 
shaft  axes.  In  the  present  case,  with  control  phasing,  if  we  represent 
the  control  system  inputs  by  OLD  A^g  and  OLD  Bjg  then  the  effective 

control  inputs  with  respect  to  shaft  axes  would  be, 

A1S'  -  OLD  A1S  cos  (<  -  £)  -  OLD  B1s  sin  (<  -  C) 

B  •  =  OLD  B1S  cos  (k  ~  O  +  OLD  A1s  sin  (<  -  £) 

If  wind  axe3  are  displaced  azimuthally  by  'I'  then  the  control 
inputs  with  respect  to  wind  axes  are. 


A1W* 

■  Ais’ 

c°3fWH  - 

B1S* 

sihTWH 

Brw' 

-  v 

C0STWH  - 

Ais’ 

sin^H 

In  the  program,  the  above  relations  are  combined  in  a  simplified 
expression  by  defining  an  effective  azimuthal  angle. 


WEFF 

=>  ? 

WH 

•f 

<  -  C 

tnus , 

Aiw’ 

-  OLD 

A1S 

COSfWEFF 

"  0LD  B,s 

Sln?WEFF 

Biw’ 

=■  OLD 

B1S 

C0SfWEFF 

+  0LD  A1S 

Sln,WEFF 

where 

tne  angles 

Aiw'  > 

0iw’ 

represent  tne 

orientation  of  tne  control 

with 

respect  to 

the  wind  axis. 

Because 

of  pitch 

flap  coupling  the 

ANF  i.3  displaced  from  the  control  axis  and  hence  the  orientation  of  the 
ANF  to  the  wind  axis  must  be  redefined.  We  have, 


A1»  *  A1W 
B1W  "  B1* 


Vi 


Vi 


where  is  the  pitch-flap  coefficient  associated  with  the  6^  ninge. 

An  additional  effect  of  blade  pitch-flap  and  pitch-lag  coupling 
is  to  modify  the  collective  pitch  input  from  the  control  system.  We  have. 


0„  -  OLD  -  ka  -  kK  £ 
c  c  3  o  5 


where  aQ,  £  are  the  steady  state  blade  flap  and  lag  angles,  and  k^  is 
the  pitch-lag  coefficient  associated  with  the  a1  hinge. 

9.  DISCUSSION 

At  thi3  stage,  only  limited  comparisons  of  model  prediction 
against  flight  trials  data  have  been  attempted.  In  the  case  of  dynamic 
response  tests  these  have  Deen  limited  to  autostabilizer  'OFF'  conditions 
where  the  aircraft  responses  are  dependent  only  on  the  aerodynamics,  and 
control  system  corrections  are  absent.  A  full  validation  program  is  under 
way  which  will  also  cover  autostabilizer  'ON'  cases  and  ASW  operations. 

Some  typical  examples  of  the  model  capabilities  and  pre¬ 
validation  limitations  are  given  in  the  following  Figures.  In  the  first 
case,  flight  data  [18]  for  steady  trimmed  flight  are  compared  with  model 
predictions.  In  Figs.  8a,b,c,  the  torque,  roll  attitude  and  pitcn 
attitude  variations  with  airspeed  show  reasonable  agreement. 

Examples  of  the  engine  model  are  given  in  Fig  9a,  b,  where  the 
transient  torque  and  rotor  r.p.m.  variations  following  a  collective  pitch 
pulse  input  are  modelled  quite  well  in  Fig  9a.  On  tne  other  nand  Fig.  9b 


(30) 


shows  that  for  a  pedal  step  input  the  engine  transients  are  much  smaller 
than  predicted  by  the  model,  although  the  yaw  rate  response  shows  good 
agreement . 

Of  the  final  examples,  Fig  10a  shows  the  variation  of  normal 
acceleration  following  a  step  change  in  collective  pitch.  In  this 
instance  much  higher  initial  loads  are  present  in  the  flight  data  because 
of  tne  time  required  for  the  rotor  inflow  to  readjust;  so-called  "inflow 
damping".  As  can  be  seen,  the  quasi-steady  rotor  represention  does  not 
cope  with  rapid  transients  during  the  first  second.  The  use  of  a  'lagged' 
inflow  in  the  model  offers  only  partial  improvement  in  predicting  the 
initial  high  loading. 

On  the  other  hand,  Fig.  10(b)  shows  results  for  a  slower  input 
(pulse).  Model  predictions  are  good  for  the  initial  portion  of  the  pulse 
but  the  need  for  the  lagged  inflow  is  evident  for  the  more  rapid 
completion  of  the  pulse.  Further  Investigations  of  the  effects  of  inflow 
lag  are  being  undertaken. 
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APPENDIX  1 


AXES  SYSTEMS 


All  axes  systems  are  right  handed  orthogonal 
Earth  (.Inertial  axes) 

The  origin  Is  any  arbitrarily  chosen  point  fixed  in  the  Earth. 

The  X£  axis  is  horizontal  along  the  intended  heading. 

The  Yg  axis  is  horizontal,  normal  to  the  X£  axis  and  positive  to 
the  right  when  looking  along  the  Xg  axis. 

The  Zg  axis  is  vertical  and  positive  downwards. 

If  desired,  the  Xg,  Yg  axes  may  be  aligned  North  and  East 
respectively. 


Body  axes 


The  origin  is  at  the  helicopter  c.g.  position. 

The  XH  axis  lies  in  the  fuselage  plane  of  symmetry,  parallel  to 
the  floor  datum,  positive  forward. 

The  Yh  axis  is  normal  to  the  Xg  axes  and  positive  to  starboard. 


The  Zu  axis  is  mutually  perpendicular  positive  downwards. 


APPENDIX  1  (CONT.) 


Shaft  axes 

As  for  body  axe3,  except  that  the  system  is  rotated  6CLI  about 

£>n 

the  body  YH  axis  so  that  0ZS  coincides  with  the  shaft  direction 
and  is  positive  downwards. 

Wind/Shaft  axes 

The  origin  is  at  the  intersection  of  the  rotor  flapping  hinge 
plane  and  the  shaft. 

The  shaft  axes  are  rotated  about  QZy  so  that  the  apparent  wind 
vector  is  contained  in  the  Xw  Zw  plane. 

The  ZH  axis  lies  along  the  shaft  in  the  direction  and  sense  of 
ZS* 

In  this  system  of  axes  blade  pitch  variation  is  given  by  the 
relationship 

ew  (Y,r)  =  e0(r)  *  A1W  COSY  -  B1w  sinY 

ANF  axes 

If  A-|W  and  B1W  are  the  components  of  cyclic  pitch  variation  with 
respect  to  the  wind  axes,  then  the  ANF  orientation  is  arrived  at 
by  successive  rotations  of  Ajy  and  -  B1W  about  the  0XW  and  the 
new  OYANF  respectively. 

Tail  Rotor  Axes  Systems 

(a)  Hub  centred  shaft  axes, 

The  XT  axis  is  parallel  to  0XU ,  positive  forward. 
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APPENDIX  2 


MATRIX  TRANSFORMS  AND  INVERSES 


The  following  rotation  matrices  and  their  inverses  are  employed 
in  the  KINEMATICS  section.  The  inverses  are  the  same  as  the  transposes. 
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FIG.  5  VARIATION  OF  YAW  RATE  STABILITY  DERIVATIVE 
(N. )  WITH  AIRSPEED. 
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■A  A»v«-Vi 


PVvXW 


V, 


vV.vV. 


DISTRIBUTION 


AUSTRALIA 

Department  of  Defence 

Defence  Central 

Chief  Defence  Scientist 

Deputy  Chief  Defence  Scientist  (shared  copy) 

Superintendent,  Science  and  Program  Administration  (shared  copy) 
Controller,  External  Relations,  Projects  and 
Analytical  Studies  (shared  copy) 

Counsellor,  Defence  Science  (London)  (Doc  Data  sheet  only) 
Counsellor,  Defence  Science  (USA)  (Doc  Data  Sheet  Only) 

Defence  Central  Library 

Document  Exchange  Centre,  DIS3  (18  copies) 

Joint  Intelligence  Organisation 
Librarian  H  Block,  Victoria  Barracks,  Melbourne 
Director  General  -  Army  Development  (NSO)  (^  copies) 

Aeronautical  Research  Laboratories 
Director 
Library 

Superintendent  -  Aerodynamics 
Divisional  File  -  Aerodynamics 
Authors:  M  J  Williams 
A  M  Arney 

R  A  Feik 
N  E  Gilbert 
C  R  Guy 
K  R  Reddy 
N  Matheson 
R  H  Perrin 
D  C  Collis 
D  A  Secomb 

Materials  Research  Laboratories 
Director /Library 

Defence  Research  Centre 
Library 

RAN  Research  Laboratory 
Library 

Navy  Office 

Navy  Scientific  Adviser 

RAN  Aircraft  Maintenance  and  Flight  Trials  Unit 
Directorate  of  Naval  Aircraft  Engineering 
Directorate  of  Naval  Aviation  Policy 
Superintendent,  Aircraft  Maintenance  and  Repair 
QIC,  Sea  King  Simulator,  RANAS  Nowra  (2  copies) 


**  PJ1  kPl-F  .-«  UK  I’*".--'  '.-v  W  ^  'A"A'^VA1STVA"AT^,.K',A*.  ^.-k*  r-n.' 


Army  Office 

Scientific  Adviser  -  Army 
Director  of  Aviation  -  Army 
Royal  Military  College  Library 

Air  Force  Office 

Air  Force  Scientific  Adviser 
Aircraft  Research  and  Development  Unit 
Library 

Technical  Division  Library 
Director  Aircraft  Engineering  -  Air  Force 
Director  of  Operational  Requirements  B  -  Air  Force 
HQ  Support  Command  (SLENGO) 

RAAF  College,  Point  Cook 

Government  Aircraft  Factories 
Manager 
Library 

Department  of  Aviation 
Library 

Statutory  and  State  Authorities  and  Industry 
Commonwealth  Aircraft  Corporation,  Library 
Hawker  de  Havilland  Aust.  Pty  Ltd,  Bankstown,  Library 

Universities  and  Colleges 
Adelaide 

Barr  Smith  Library 

Flinders 

Library 

La  Trobe 
Library 

Mel  bourne 

Engineering  Library 
Monash 

Hargrave  Library 

Newcastle 

Library 

Sydney 

Engineering  Library 
Professor  G  A  Bird 
Mr  J  Blackler 


NSW 

Physical  Sciences  Library 

Professor  R  A  A  Bryant,  Mechanical  Engineering 
Professor  G  D  Sergeant,  Fuel  Technology 


Queensland 

Library 

Tasmania 

Engineering  Library 

Western  Australia 
Library 


RMIT 

Library 

NRC 

Aeronautical  &  Mechanical  Engineering  Library 


France 

ONERA,  Library 
India 

Hindustan  Aeronautics  Ltd,  Library 

National  Aeronautical  Laboratory,  Information  Centre 

United  Kingdom 


Royal  Aircraft  Establishment 
Bedford,  Library 

British  Library,  Document  Supply  Centre 
Westland  Helicopters,  Limited 

United  States  of  America 

NASA  Scientific  and  Technical  Information  Facility 


SPARES  (15  copies) 
TOTAL  (105  copies) 


Department  of  Defence 


DOCUMENT  CONTROL  DATA 


U.ARNa  I  3>.  Establishment  No 

AR-004-467  ARL- AERO- TM- 379 

4.  Title 

A  MATHEMATICAL  MODEL  OF  THE 
SEA  KING  MK.50  HELICOPTER 
AERODYNAMICS  AND  KINEMATICS 

8.  Authors) 

M.J.  WILLIAMS 
A.M.  ARNEY 


10.  Corporate  Author  anil  Address 


Aeronautical  Research  Laboratories 
P.0.  Box  4331 , 

MELBOURNE,  VIC.  3001 


12.  Secondary  Distribution  (of  this  document) 

Approved  for  public  release. 


2 


.  Document  Dote 


3.  Todt  No 


JUNE  1986 

5.  Security 
a.  document 

UNCLASS 


DST  85/028 

6.  hbPoge* 

44 


b.  title  c.  abstract 


7.  No  Kefs 


_U 

9. 


JJ  _  _ 

Downgrading  Instructions 


19 


1 1.  Authority  (os  appropriate) 

□.Sponsor  b. Security  c .Downgrading  d .Approval 


Overseas  enquirers  outside  stated  limitations  should  be  referred  through  ASDIS,  Defence  Information  Services  Branch,  Department  of 
Defence,  Campbell  Pork,  CANBERRA  ACT  2601 

13d.  This  document  may  be  ANNOUNCED  in  catalogues  and  awareness  services  available  to  — 

No  limitations  . 


1 33).  Citation  for  other  purposes  (ie  casual  announcement)  may  be  (sete  :l)  unrestricted(or)  as  for  13  o. 
14.  Descriptions 

Aerodynamics  Mathematical  models 

Kinematics  Helicopters 

Sea  King 


IS.  COS  ATI  Group 
01030 
01010 


76.  Abstract 


> Details  are  given  of  the  expressions  used  to  describe  the 
aerodynamics  and  kinematics  of  the  Sea  King  Mk.50  helicopter  during  steady 
flight  and  low  rate  manoeuvres  up  to  an  advance  ratio  of  0.3.  The 
aerodynamics/kinematics  formulation  is  a  major  component  of  the  Sea  King 
mathematical  model  developed  by  Aeronautical  Research  Laboratories  (ARL)  for 
flight  simulation  of  this  Anti-Submarine  Warfare  helicopter. 
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